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Dynamics of central collisions of heavy nuclei in the energy range from few tens 
of MeV/nucleon to a couple of GeV/nucleon is discussed. As the beam energy 
increases and/or the impact parameter decreases, the maximum compression in- 
creases. It is argued that the hydrodynamic behavior of matter sets in the vicinity 
of balance energy. At higher energies shock fronts are observed to form within 
head-on reaction simulations, perpendicular to beam axis and separating hot com- 
pressed matter from cold. In the semicentral reactions a weak tangential disconti- 
nuity develops in-between these fronts. The hot compressed matter exposed to the 
vacuum in directions parallel to the shock fronts begins to expand collectively into 
these directions. The expansion affects particle angular distributions and mean 
energy components and further shapes of spectra and mean energies of particles 
emitted into any one direction. The variation of slopes and the relative yields mea- 
sured within the FOPI collaboration are in a general agreement with the results 
of simulations. As to the FOPI data on stopping, they are consistent with the 
preference for transverse over the longitudinal motion in the head-on Au -|- Au 
collisions. Unfortunately, though, the data cannot be used to decide directly on 
that preference due to acceptance cuts. Tied to the spatial and temporal changes 
in the reactions are changes in the entropy per nucleon. 

1 Introduction 

In my talk, I shall discuss the dynamics of central heavy-ion reactions in the 
energy range from few tens of MeV to a couple of GeV per nucleon. I shall use 
simple estimates, data, and transport reaction simulations, to reach conclu- 
sions. When moving from low to high energies, from light to heavy systems, 
or from high to low impact parameters, I will argue, the onset of hydrody- 
namic behavior in collisions is demarked by the balance energy. Transport 
simulations indicate that in the high-energy head-on collisions of nuclei, in the 
hydrodynamic regime, nuclear shock waves develop, perpendicular to the beam 
axis, separating hot compressed nuclear matter from cold matter. With time, 
the compressed matter exposed to the vacuum acquires collective motion giving 
rise to a characteristic behavior of transverse cm spectra of different particles 
with mass. At finite impact parameters a so-called weak discontinuity develops 
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in-between the shocks. Both the motion associated with the discontinuity and 
the expansion contribute to the sideward flow. Recently FOPI data became 
available on stopping and on expansion in central collisions. I shall examine to 
what extent these data are consistent with the above scenario. Finally, I shall 
discuss the production of entropy in collisions and conclude. 

2 Onset of Hydrodynamic Behavior: Balance Energy 

The hydrodynamic behavior requires a local equilibrium and, thus, a small 
mean free path, A ^ l/(/9crNN) ~ 2 fm, compared to the size of a system, 
X/R <C 1. Here p is density and the effects of Pauli principle are ignored. 
The condition implies that the hydrodynamic behavior is more likely in the 
interactions of heavy than of light nuclei, and in the more central than in 
peripheral collisions. Adding to this consideration, the effects of the Pauli, 
lengthening the mean free path, can slow down the equilibration required for 
the hydrodynamic behavior, at low beam energies. At moderate energies the 
momentum dependence of nuclear mean field may cause a reduction in the 
number of NN collisions affecting equilibration. 

At low energies, nuclei colliding at moderate to high impact parameters 
are expected to clutch and rotate in a binary system, emitting particles to 
negative angles. At high energies, high pressure should develop in the region 
between the nuclei, expelling nucleons and the residual nuclei to positive angles. 
For emission quantified in terms of the average momentum in the reaction plane 
as a function of rapidity the slope of the dependence on rapidity 

should then change sign, when going from low to high bombarding energies, 
at the so-called balance energy. I will argue that, whether in the bombarding 
energy or in the impact parameter, the balance energy demarks the onset of 
the hydrodynamic behavior. 

By investigating the magnitude of the flow, d{p^/A)/dy at midrapidity, 
the balance energy has been studied El as a function of the reduced impact 
parameter h — b/bmax, in the system ^"Ar -I- ^^Sc. The energy was observed 
to rise with the rise in b, see Fig. |l|. The measurement is consistent with the 
hydrodynamic behavior setting earlier in the more central than in the more 
peripheral collisions as the beam energy increases (see also Rcf . u) . The depen- 
dence of the balance energy on total mass for symmetric systems is shown in 
Fig. 1^. The balance energy is higher for light systems and lower for heavy sys- 
tems, consistent with the hydrodynamic behavior setting earlier in the heavy 
than in the light systems. The balance energy calculated with the transport 
modelB strongly depends on NN cross sections; better agreement with the data 
is obtained for reduced than for free-space cross-sections. 
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In the measurements such as e.g-PflB there has been no control on the di- 
rection of the reaction plane. Thus only a suppression of the slope of {p^ /A){y) 
at midrapidity could be observed as a function of the beam energy. Principally, 
one might imagine that at the deduced balance energy only a suppression of 
the slope takes place without a change in sign. Recently in the experiment^ 
though, the polarization of 7 rays from a residual target nucleus has been mea- 
sured, in the ^"^N -I- ^^''Sm reaction at 35 and 100 MeV/nucleon, in coincidence 
with the emission of energetic a particles at 9iab — 30°. For intermediate 
reduced impact parameters 6, the 7 polarization was found to change sign in- 
between 35 and 100 MeV/nucleon, in a manner consistent with the change of 
sign of the flow, from negative (attractive) to positive (repulsive). For high b, 
the sign of the polarization was found such as for an attractive flow at either 
of the beam energies. 

Let us try to understand the mass dependence of the balance energy in 
some detail; we shall consider a symmetric reacting system. As the nuclei come 
into contact, the forces that act on each include the force due to pressure in the 
neck region with equilibrium and dissipative terms, the proximity force acting 
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Figure 2: Balance energy in nearly symmetric systems as a fumjtion of total mass. Filled 
squares and circles represent, respectively, data of Ref. □ and □ Open symbols represent 
transport-model calpiilations of Ref. □ The dashed and dotted lines indicate the general 
trends of the resultsu The solid line represents an estimate according to Eq. (O). 



within the circumference of the neck, and the Coulomb force, 

F = Foq+dis + Fprox + Fcoul • (1) 

The pressure within the neck region is 

P = Po^ + ^V^^~2dp\/^(^£^ , (2) 

where the first term on the r.h.s. is equilibrium pressure, next is a portion of 
the viscosity tensor that contributes to transverse force, and final term is the 
finite-range correction. The factor rj in the viscous term is viscosity coefficient. 
The proximity force may be estimated either from the finite-range term in 
or by considering the surface energy within the system, obtaining 

-Fprox = -2^ i? a, (3) 

where cr ~ 1.0 MeV/fm^ is the surface tension and R is the nuclear radius. 
The proximity force is proportional to neck circumference. The Coulomb force 
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may be estimated as Fcoui = ZiZ2eV(i?i + i?2)^ ^ 0.016 A^/^ MeV/fm. 
Given that no compression is expected at high impact parameters^ the equi- 
hbrium pressure can be estimated as Pcq = 2po ^'cm = Po-£'iab/6. Finally, 
the dissipative viscous correction to pressure may be estimated as Pdis = 
— (2/3) {poPf '^/S) vuh/R (an additional, nominally dissipative, correction may 
arise from the momentum dependence of the mean field) . At the balance en- 
ergy the net force acting on any of the nuclei vanishes, i.e. 

= Fprox + Fcoul + TT i?2 (P,q + Pdis) , (4) 

where tt has been used as an estimate for the neck area. Solution of Eq. (||) 
gives 

E^^, 7e^ f4 V + ^ - 0.21 A^/^ 



^RJ Rpo 
740 160 



^2/3 ^1/3 



0.21^2/3, (5) 



where Ci;- is Fermi energy and the last result is given in MeV {A is the total mass 
of the system). The crude estimate is shown with a solid line in Fig. (H). 
Quantitative result allows to conclude the following. The nonequilibrium cor- 
rection in the pressure is very important; if this correction were omitted the 
balance energy would turn out to be 4 times too low compared to data. Ob- 
served large balance energies demonstrate that little compression is reached 
in periphery; if one assumed that all cm. energy per nucleon were used up 
in a compression, the balance energy from (^) would underestimate data by 
an order of magnitude. The sensitivity to elementary cross sections for light 
nuclei, Ebai oc a^'^ with a^n < o^np, seems to explain the rise of balance energii. 
with the rise in the isospin asymmetry found in the transport calculations; □ 
the weakening of the mean field should act to reduce rather than to increase 
the balance energy. 



3 Head-On Reaction Simulations 

Dynamics of nuclear reactions above the balance energy has been studied in 
some detail following simulations based on the Boltzmann equation with a den- 
sity and isospin dependent mean field, free NN cross-sections, and composite 
production, in Ref. u and elsewhere. It is possible to break the dynamics 
into several stages. A 400 MeV/nucleon Au + Au system at 6 = 0, for which 
baryon-density contour-plots are shown in Fig. ^, will serve to illustrate points. 
Unless otherwise indicated the results of simulations will correspond to a stiff 
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Figure 3: Contour plots of baryon density in the reaction plane in Au + Au collisions at 
400 MeV/nucleon. The displayed contour lines are for the densities p/pa = 0.1, 0.5, 1, 1.5, 
and 2. 



equation of state with an incompressibility oi K — 380 MeV. Following an ini- 
tial interpenetration of projectile and target densities, the NN collisions begin 
to thermalize matter in the overlap region making the momentum distribution 
there centered at zero momentum in the c.m.s. The density in the overlap 
region rises above normal and a disk of excited and compressed matter forms 
at the center of the system. More and more matter dives into the region with 
compressed matter that begins to expand in transverse directions. At late 
stages, when the whole matter is excited, transverse expansion predominates. 
A further view at the situation in reaction at 6 = is given in Fig. ^. 

If a nuclear system were very large, then it would be necessarily ruled by 
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Figure 4: Baryon density (top panels), radial velocity (center panels), and entropy per 
baryon along (solid lines) and perpendicular (dashed lines) to the beam axis at different 
indicated times, in the 6 = collision at a beam energy of 400 MeV/nucleon. 



laws of hydrodynamics. In the initial state of a head-on reaction, a disconti- 
nuity exists in the velocity field at the contact area of two nuclei, with the ve- 
locities in the nuclei directed opposite. Within hydrodynamics, such an initial 
discontinuity would have to break at finite times into two shock fronts travel- 
ling in opposite directions into the projectile and target. To assess whether the 
interfaces between normal and compressed matter, cf. Figs. ^ and Q may be 
interpreted as the shock fronts, one can ask whether the state of the matter at 
the center is such as expected behind a developed front. Properties of the mat- 
ter behind a front can be determined, in a straightforward manner, from the 
conservation of hydrodynamic fluxes of baryon number, momentum, and en- 
ergy. Figure |^ displays some of these properties, baryon density, temperature, 
and entropy per baryon, together with a maximum density and entropy per 
baryon at a maximum density, from simulations. The latter quantities gener- 
ally follow the trends anticipated from hydrodynamics. The shocks are nearly 
completely developed in the high-energy simulations of the head-on reactions 
of heavy nuclei. 
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Figure 5: Solid lines show the baryon density (top panel), temperature (center panel), and 
entropy per baryon (bottom panel) expected behind a developed shock front at 6 = 0, as a 
function of beam energy, from solving the Rankine-Hugoniot equation. Filled circles indicate 
the maximum density (top panel) and entropy per baryon at maximum density (bottom 
panel) from the simulations of fe = Au + Au collisions at different beam energies. Crosses 
in the top panel indicate the maximum density from the simulations of 1 GeV/nucleon Au 
+ Au collisions at b = 3, 6, 8, 9, 10, 11, and 12 fm, plotted against the beam energy scaled 
down by a factor equal nonrelativistically to (1 — 6-^/4/?^), cf. Eq. (H). 

Let us now turn attention to the expansion. The hot matter in-between 
the shocks is exposed in transverse directions to vacuum and, as the pressure 
of the hot matter is finite and that of vacuum is zero, the matter begins to 
expand coUectively into transverse directions. The features of this process may 
be understood at a qualitative level in terms of the self-similar cylindrically- 
symmetric hydrodynamic expansion. In the self-similar expansion, the velocity 
is proportional to the distance from symmetry axis, v = T{t) r, which roughly 
holds in the direction of 90° in the Au -|- Au collision at 400 MeV/nucleon, 
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at the center of the system, see Fig. ^ On solving, in a near-analytic manner, 
the hydrodynamic equations of motion for the expansion^ one arrives at the 
expectation that the matter should accelerate collectively in the transverse 
direction for a time R/cs, where Cs is the speed of sound. This time is such as 
necessary for a signal to propagate into the interior. It is generally consistent, 
R/Cg ~ 20 fm/c, with the time during which collective transverse energy, 
calculated from the local collective velocities 

Ei,, = J drdp/mjv(«^)V2, (6) 

rises in the simulation displayed in Figs. || and ^. 

After shocks reach the vacuum after a time ^ {2R/vo){pi — po)/pi ~ 
15 fm/c, where vq is the initial velocity in the cm. and pi is the density 
in the shocked region, an expansion along the beam axis sets in, see Fig. ^ 
However, as the expansion in transverse directions is already in progress and 
matter becomes decompressed, the expansion in longitudinal direction does 
not acquire same strength. 

While definite features of reaction dynamics appear to be consistent with 
a hydrodynamic behavior of the matter, there are also important differences. 
For example, in a system continuing to expand hydrodynamically temperatures 
would have dropped to zero and all kinetic energy would have got converted 
into collective energy of expansion. For that collisions would need to continue 
down to very low densities. The transverse collective energy from (^) freezes 
out at about 40 MeV/nucleon. With about 10 MeV/nucleon of longitudinal 
collective energy, one finds that for nucleons the net collective energy con- 
stitutes about 50% of the kinetic energy in the cm. for a 400 MeV/nucleon 
reaction. 



4 Semicentral Reactions 

In the initial state, an angle a of inclination of the plane of discontinuity in 
velocity, relative to the beam axis, is given nonrelativistically by 

cosaw^. (7) 

While at 6 = velocities are normal to the plane of discontinuity, the veloci- 
ties have finite tangential components relative to that plane at finite b. Such 
components are continuous across the shocks which detach from the initial 
plane of discontinuity and, with tangential components being directed oppo- 
site for the matter from a projectile and from a target, at the center of the 
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Figure 6: The initial discontinuity between the projectile and target velocities breaks at 
a finite b into two shock fronts propagating into the projectile P and target T (thick solid 
lines), and a weak tangential discontinuity in-between (dotted line). Shock- front position 
within projectile at a later time is indicated with a dashed line. 

system a so-called tangential discontinuity develops, see the sketch in Fig. |[ 
Density, pressure and entropy are continuous across that discontinuity. Only 
the tangential velocity component changes. 

As only the normal velocity component drops to zero across shocks, the 
shocks are weaker at a finite b than in a head-on reaction. With the normal 
component being equal to vo sin a, the effective cm. kinetic energy for a shock 
at a finite b becomes 

_ . sm « = — (1 - COS a) = — (^1 - — j , (8) 

rather than Eiab/4: in the nonrelativistic approximation. In consequence, the 
density behind a shock at a finite b should be lower than at 6 = and it should, 
in fact, coincide with the density in a & = reaction at beam energy reduced 
by a factor (1 — b^/AR^). That is tested in the upper panel of Fig. |[ In the 
1 GeV/nuclcon reaction the density appears to follow the above expectation 
up to very high impact parameters. 

While a shock width generally does not change with time, if conditions at 
the two sides of a shock wave stay the same, the width of a weak discontinuity 
increases with time. That this is the case, may be seen by adapting the Navier- 
Stokes equations to the conditions in the vicinity of a weak discontinuity and 
by observing that the equations reduce there to a diffusion equation for a tan- 
gential velocity component. The diffusion coefficient is equal to the kinetic 
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Figure 7: Baryon density p (top panel), normal and tangential velocity-components, ti„ and 
vt, to the discontinuities in nuclear-matter (center panel), and entropy per baryon (bottom 
panel) as a function of the distance from the center of a 400 MeV/nucleon Au + Au system 
at 6 = 7 fm, along the normal to the discontinuities. At the given time t = 13.5 fm/c the 
discontinuities (at their centers) are inclined at an angle a ~ 52° relative to the beam axis. 



viscosity coefficient v = ri/{miqp)^ where rj is viscosity. With an initial condi- 
tion for the tangential velocity in the ideal-fluid limit in the form vt = 
the solution to the diffusion equation becomes vt {r, t) — erf (r/2\/t4) , where 
r is the distance along a line perpendicular to the discontinuity, i.e. the width 
increases as t^/"^. If a nuclear system were very large, then the two shock waves 
and a weak discontinuity would well separate from each other, as the distance 
between the shocks increases as t. In actual simulations, though, the tails of 
discontinuities overlap, cf. Fig. ^ In that figure, the center of the shock is at 
a position r ~ 3.3 fm, where density rises to the mean of that in normal matter 
and that at the center of the system, where, further, entropy per nucleon rises 
to a half of its maximum value, and where magnitude of the normal velocity 
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component drops to a half of the asymptotic value. The center of the tangen- 
tial discontinuity in the figure is, on the other hand, at r = 0, where tangential 
velocity component vanishes. The width of the discontinuity is characterized 
by (5 1.7 fm which agrees well with S = ^JrftJniNp expected from the sim- 
ple formula, given above, for the space and time behavior of the tangential 
velocity. 

Let us start to confront the findings with some data. Just as at 6 = 
in a heavy system, an expansion develops at a finite h in-between the shock 
waves, as an effect of the exposure of the hot matter to vacuum. Anisotropy 
in the expansion, with regard to the direction perpendicular to the reaction 
plane and the direction of the shock motion, results from the delay in the start 
of expansion in the latter direction, as in the case 6 = 0. The anisotropy 
(squeeze-out) may be quantified with a ratio of the eigenvalue of the kinetic 
energy tensor associated with the direction out of the reaction plane, to the 
smaller of the eigenvalues associated with the directions within the plane, 

^ If. (9) 

As 6 —> in a heavy system (but not in light), the ratio R21 tends to {E^) /2(£'ll). 
The Au -I- Au data for the ratioEl are compared to the results of calculations, 
corrected for the experimental inefficiencies, in Fig. |8| At 250 MeV/nucleon 
and 400 MeV/nucleon a degree of agreement between the data and the calcu- 
lations is observed. The inefficiencies reduce the sensitivity of the ratio to the 
nuclear compressibility. 

With regard to the directions within the reaction plane, both the tangential 
motion behind the shock waves for finite impact parameters, and the anisotropy 
in expansion, contribute to the sideward deflection quantified usually in terms 
of the slope of mean momentum within the reaction plane at midrapidity, 
F = d{p^ /m)/dy. Measured and calculated slopes in the 400 MeV/nucleon 
Au 4- Au reaction are shown in Fig. ||. If the motion associated with the weak 
discontinuity were not contributing to the sideward deflection, then the slope 
would have been limited from aboveQ by x/(2-\/l + x) where x = ((-^^)/2 — 
(ijll))/ (£'11 ) and the energy components are from a 6 = reaction. For protons 
from 400 MeV/nucleon reactions this limit is equal to 0.22, while slopes given 
in Fig. ^ reach values twice as high. 

5 Comparisons to FOPI Data 

As 6 — > 0, the expansion of matter between the shocks, favoring the transverse 
direction of motion, competes with the transparency effects in the corona, 
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Figure 8: Ratio of the out-of-reaction-plane mean energy-component to the lower in-plane 
eigenvalue of the kinetic energy tensor, as a function of reduced impact parameter, in 
Au + Au coUisioHfi at beam energies indicated in the panels. Circles represent the mea- 
surements of Ref. □. Solid and dashed lines represent the results of the transport-model cal- 
culations with mean fields yielding a stiff and soft equation of state, respectively. The dotted 
lines represent the results of the cascade model. 



that favor the longitudinal direction. For light systems in simulations, such 
as Ca -I- Ca, the transparency effects prevail, leading to (E-^) /2{E^ < 1 at 
6 = 0, but not so for heavy systems such as Au -I- Au, as was discussed and 
as is evidenced in Fig. ^. The experimental detection of (E^) /2{E^ > 1, 
i.e. (E-^) / {E^ > 2, would demonstrate a violence of the nuclear hydrody- 
namic phenomena in heavy systems. On the other hand, if for any reason 
the in-medium NN cross sections were lower than the free-space cross-sections 
used in the simulations, then the transparency effects could prevail over the 
anisotropy in the expansion, even within heavy systems. A clear experimental 
answer requires a full 2t: coverage either in the forward or in the backward 
cm. hemisphere. 

Generally, the experimental determination of the average anisotropy of 
momentum distribution in 6 ~ collisions may be difficult, due to the contam- 
inations of event samples by events cacresponding to intermediate b. The idea 
borne within the FOPI collaboration til has been then to employ, for compar- 
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Figure 9: Left panel shows the flow parametp^ F in Au + Au reactions at 250 and 400 
MeV/nucIeons, from the measurements of Ref. tj, as a function of the reduced participant 
proton multiplicity. Right panel shows the flow parameter from the simulations of Au + Au 
reactions at 400 MeV/nucleon, as a function of impact parameter. 



isons to model calculations, the differential cross sections for anisotropics quan- 
tified in terms of the event-by-event energy ratio, sometimes termed ERAT, 

(10) 

where the sums extend over detected charged particles from collisions emerging 
in the forward cm. hemisphere. In Phase I, the FOPI detector setup covers 
laboratory angles 1° < diab < 30° for the full azimuth, with a hroa.ijer range 
expected in Phase II. The cross sections for the energy ratio measuredE5E3ild in 
the 250 MeV/nucleon Au + Au reaction are shown in Fig. |l^ together with the 
results of calculations filtered through the acceptance of the FOPI apparatus. 
The calculations were carried out for a stiff (H) and for a soft (S) equation 
of state, without {H and S) and with (SM) momentum dependence, for free- 
space cross sections and for cross sections reduced (in the H case) by a factor 
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Figure 10: Differential ERAT (cf. Eq. (0)) cross scctiorii in Aii + Au reaction at 
250 MeV/nucleon as a function of ERAT. The data of Refs. BliaM are represented by 
the histogram. The results of calculations utilizing a momentum-independent mean field, 
that yields a stiff equation of state, are represented by squares for free NN cross sections, and 
by stars for density-dependent NN cross sections with v = 0.15, cf. Eq. (|ll[). The results of 
calculations utilizing free cross sections and mean fields, that yield a soft equation of state, 
are represented by crosses and diamonds, respectively, for fields with and without momentum 
dependence. Final momentum distributions in the calculations have been filtered through 
the experimental acceptance. 

exponentially dependent on density, 

(Jnn = ^NN exp(-i^p/po) • (11) 

Because of the acceptance cuts, no events above or in the vicinity of 
E±/E^^ — 2 can be actually observed in Fig. |l^. The comparison between 
the data and calculations at 250 MeV/nucleon and at other energies, shows 
that the data allow only for a small reduction {i/ ~ 0.05 — 0.15) in the NN 
cross-sections in the simulations with a stiff equation of state and hardly for 
any reduction in the simulations with a soft equation. Values as large as > 0.4, 
for the stiff equation of state, are required for the transparency to prevail over 
the transverse expansion in the b — Au -f Au collisions. When the acceptance 
cuts are removed, the ERAT cross sections from simulations extend well above 
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Figure 11: Differential ERAT cross section in Au + Au reaction at 250 MeV/nucleon 
as a function of ERAT. The data are represented by the histogram. The results of an H 
calculation utilizing free NN cross-sections are represented by open and filled squares for 
the cases with and without experimental cuts, respectively. The symbols '0', '1', and '2', 
represent contributions to the ERAT cross section obtained without experimental cuts, 
from the impact-parameter intervals < 6 < 0.5 fm, 0.5 fm < 1.5 fm, and 1.5 fm < 2.5 fm, 
respectively. All the results of calculations were obtained using the forward cm. hemisphere 
only. 



ERAT = 2, see Fig. ^ It should be possible to observe such distributions 
within the Phase II of the FOPI experiment. Figure 12 shows the differential 
ERAT cross section computed with cuts and broken down into contributions 
from narrow 6-intervals. The figure demonstrates that selecting the highest 
ERAT values can yield a surprisingly pure sample of most central events with 
6 < 1.5 fm. 

The transport simulations demonstrate various effects of the anisotropy in 
the collective expansion on observables from the collisions. However, the collec- 
tive expansion should also affect observables at any one angle. Thus, in a glob- 
ally equilibrated system, exponential spectra would be expected at any angle, 
with the same slope parameter for different particles, equal to the tempera- 
ture, dNx/dp^ (X exp(— i?a;/r). In the presence of the collective expansion, on 
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Figure 12: DifTerential ERAT cross section in Au + Au reaction at 250 MeV/nucleon 
as a function of ERAT. The data are represented by the histogram. The results of an H 
calculation with free NN cross-sections are represented by squares. The symbols '0', '1', 
'2', . . . , 'A', represent contributions to the calculated cross section from impact parameter 
intervals < 6 < 0.5 fm, 0.5 fm < 1.5 fm, and 1.5 fm < 2.5 fm, . . . , 9.5 fm < 10.5 fm, 
respectively. 

the other hand, the particle spectra should exhibit different slopes for particles 
with different mass, with particle spectrum becoming flatter as particle mass 
increases. This is due to an increased sensitivity to the collective motion with 
an increasing mass; for the average kinetic energy of a particle, one would, 
in fact, expect a linear rise with the mass number, 

{e.) = It+!!^^It + a^!^, (12) 

where v is collective velocity and a uniform temperMure T is assumed. The spec- 
tra of light particles from the FOPI measurementsliJ and from the calculationsQ 
are shown in Fig. ^|, together with the average energies as a function of the 
mass number and with the relative particle yields. Similar flattening of the 
spectra with particle mass is observed for the data and for the calculations. 
The average kinetic energies are consistent and they rise approximately linearly 
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Figure 13: G«mparison of light-particle measurements!^ (filled symbols) to transport-model 
calculations lI (open symbols) for the central Au + Au reaction at 250 MeV/nucleon. Left, 
upper right, and lower right panels show, respectively, transverse cm. spectra, mean cm. en- 
ergies, and transverse particle yields normalized to proton yields. Thin straight line across 
the cm. spectrum is drawn to guide the eye, as parallel to the helion spectrum. 



with the mass for the hghtest fragments, exhibiting the presence of the col- 
lective expansion. Also the particle yield ratios, shown in Fig. are rather 
consistent at 250 MeV/nucleon between the theory and experiment. These 
ratios reflect the entropy produced within the system, as, in the limit of Boltz- 
mann statistics, 

^ cx exp ^-Aj. . (13) 
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Figure 14: Left panel shows entropy per nucleon as a function of reduced participant 
multmlicity, determined from 400 MeV/nucleon Au + Au wide-anrt light-particle data of 
Ref . lij (stars) and from intermediate- mass fragment-idata of Ref. E3 (cross). Right panel 
shows, as a function of impact parameter, calculated Q total entropy per nucleon associated 
with nucleons emitted into wide angles (circles and dashed line) and contributions to that 
entropy from shock waves (solid line) , from diffusion of vorticity around weak discontinuity 
(between the solid and dashed lines), and from expansion of participant region (between the 
dashed and dotted lines). 



6 Entropy Production 

While the entropy produced in reactions attracted much experimental atten- 
tion in the past, this has not been the case recently. Yet, the entropy ties in 
an interesting fashion to the reaction dynamics and it may, principally, provide 
information on the low-density phase transition. 

The entropy in central reactions of heavy nuclei comes from several sources. 
Most of the entropy per nucleon in the participant region at 6 ~ is produced 
in shock waves, see Figs. ^, ||, and |lj. Additional entropy is produced during 
expansion. The hydrodynamic Navier-Stokes equations yield for the rate of 
change of the local entropy per baryon with time 



d_ 

di 



V 
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1 , . r] f dv, dvk 2 dvA 
= -d:v (.VT) + 2^ (^^ + ^ - 3 ^ j , (14) 

where k is heat capacity coefficient and 77 is viscosity coefficient. Whereas 
the viscosity term on the rhs of always leads to the rise in local entropy, 
the heat conduction term may cause a local drop, although still the entropy 
increases on a global scale. The latter is the case when the temperatures on 
the average decrease as one moves away from a considered point. 

In the collective expansion, the edges of participant matter cool off; in 
an isentropic expansion pT^^I^ = const and thus T as p — > 0. Then 
the heat conduction, or simply flux of energetic nucleons from the inner re- 
gion, reheats the outside matter, significantly raising the entropy there, see 
Fig. I Within the inner region the entropy per nucleon decreases. For & ~ 
400 MeV/nucleon Au + Au collision, analytic estimates with Eq. ( p^ yield 

an overall increase in the entropy 5{S/A) ~ 2-^ ( pj"'"'' J ~ 1 on account of 
the heat conduction, and, practically, no increase on account of viscosity due 
to the lack of shear in a central collision. As b increases, shear appears and 
then the entropy gets to be produced through viscosity as the weak viscos- 



ity spreads, cf. Fig. 14, Moreover, the production of the entropy per nucleon 
in the participant region during expansion rises, since the size of the region 
diminishes, 5{SIA) oc A'^/'^.^ 

In a recent investigation^, Petrovici et al. attempted to describe, simul- 
taneously, the yields and the average energies of intermediate-mass fragments 
from central 250 MeV/nucleon Au -f Au collisions, in terms of an analytic 
ideal-fluid hydrodynamic model l3 combined with a model of statistical dis- 
assembly. Within an ideal-fluid model the entropy per nucleon is conserved 
along a streamline; there is no entropy production. If the initial entropy rises 
or declines with distance the from the center, so does the entropy towards the 
freeze-out. The parameter within the model, regulating the shape of initial 
density, regulated simultaneously the entropy. Best fit to data was obtained 
for the entropy sharply rising with the distance from the center, cf. Figs . |l5| 
and[l^, consistent with Fig. ^. Such variation of the entropy favors, cf. Eq. ([l3|), 
the production of heavier fragments within the inner region of matter and of 
lighter fragments outside, see Fig. The collective energy per nucleon is lower 
inside and higher outside within the model and according to the simulations. 

As is apparent, the change in entropy with the impact parameter is tied to 
quite nontrivial aspects of reaction dynamics. Even at one impact parameter, 
in very central reactions, fragments with different masses may be associated, 
due to a variation of the entropy, with different regions and emission times. 
It may be worthwhile to study directly the entropy and the collective energy 
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Figure 15: Initial fireball density (top left panel), temperature (top right), and entropy 
per nucleon (bottom left), and final fragment multiplicity (bottom right) as a functimpi 
position within the fireball. Solid and dotted lines correspond to the shape parameter EZI'LJ 
a = 3.0 and a = 0.01, respectively. 

as a function of mass accessing the spatial and temporal changes in a reaction. 
7 Conclusions 

To summarize, the dynamics of central energetic heavy-ion collisions proceeds 
through numerous complex stages. As the beam energy increases and/or the 
impact parameter decreases, the maximum compression increases. Hydrody- 
namic behavior of matter sets in the vicinity of balance energy. At that energy, 
the bulk pressure competes with viscous effects and with the proximity force. 
The 7 polarization measurements show directly the change in sign of the nu- 
clear deflection at the balance energy. At higher energies shock fronts are ob- 
served to form in the head-on reaction simulations, perpendicular to beam axis 



21 



N 



> 

0} 



< 

LLI 

V 




5 6 T 



Figure 16: (a) Avepage multiplicity in the forward hemisphere of central 250 MeV/nucleon 
Au + Au reaction tj and (b) average kinetic energy per nucleon for different azimuthal 
intervals, as a function of fragment charge. Measurements are represented by open symbols. 
Results of the model calculation for the optimal a = 0.03 filtered through the experimental 
acceptance are represented by filled triangles. 



and separating the hot compressed matter from the cold. In the semicentral 
coUisions a weak tangential discontinuity develops in-between the shock fronts. 
The hot compressed matter exposed to the vacuum in directions parallel to the 
shock fronts begins to expand collectively into these directions. The expansion 
affects particle angular distributions and mean energy components and further 
shapes of spectra and mean energies of particles emitted into any one direction. 
The variation of slopes and the relative yields measured within the FOPI col- 
laboration are in a general agreement with the results of simulations. As to the 
FOPI data on stopping, they are consistent with the preference for transverse 
over the longitudinal motion in the head-on Au -I- Au collisions. However, the 
data cannot be used to decide directly on that preference due to acceptance 
cuts. Tied to the spatial and temporal changes in the reactions are changes 
in the entropy per nucleon. Different fragments are associated with regions of 
different entropy values. Besides the participant entropy varies with impact 
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parameter. The entropy sources include shock waves, weak discontinuity and 
dissipation in expansion. 
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